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We develop a model describing the electrochemical conversion of an organic semiconductor (specif- 
ically, the active material in a light-emitting electrochemical cell) from the undoped non-conducting 
state to the doped conducting state. The model takes into account both strongly concentration- 
dependent mobility and diffusion for the electronic charge carriers and the Nernst equation in the 
doped conducting regions. It is demonstrated that the experimentally observed doping front pro- 
gression in light-emitting electrochemical cells can be accurately described with this model. 



The pioneering demonstration that it is possible to per- 
form doping on an organic conjugated polymer and attain 
a high metallic-like electronic conductivity was awarded 
with the Nobel Prize in Chemistry in 2000 [ij . The gen- 
eral opportunity for a controlled tuning of the electronic 
and optical properties of organic semiconductors via var- 
ious doping techniques has opened up for a wide range of 
emerging novel and flexible applications. It also provides 
for important and stimulating science as regards to the 
fundamental processes in organic semiconductors. In the 
latter context, it is relevant to find and understand the 
differences and similarities to the well-established physics 
and chemistry that is dictating the behavior of conven- 
tional inorganic semiconductors. 

A striking difference between organic and inorganic 
semiconductors is that the former are soft materials with 
weak intermolecular bonds and low dielectric constants, 
and as such are strongly influenced by electron-lattice in- 
teractions [2] • The mere existence of an excited electronic 
state on an organic semiconductor creates a significant lo- 
cal lattice distortion, which causes the electronic charge 
carrier to self-localize over a spatial volume of nm-sized 
dimensions (and to be termed an electron/hole polaron). 
Moreover, local disorder is common in polymers, and the 
existence of a random distribution of dopant counter-ions 
will create even more energetic disorder and further lo- 
calize the electronic charge carriers 0, Q . 

A number of interesting features of organic semicon- 
ductors originate in their specific properties and distin- 
guish them from inorganic semiconductors. First, elec- 
trochemical (and also chemical) doping can be performed 
with straightforward means in-situ. By applying an ap- 
propriate voltage to an electrode coated with an organic- 
semiconductor film and in contact with an electrolyte, 
it is possible to inject electronic charge into the organic 
semiconductor, which subsequently is electro-statically 
compensated by injection (ejection) of ions from (into) 
the electrolyte The necessary motion of ions within 

the organic film is facilitated by its soft, and in some 
cases porous, nature. Second, the doping levels corre- 



lating to a high electronic conductivity is much higher 
in organic semiconductors than in their inorganic coun- 
terparts (~0.1 vs. ~10~'' dopants/repeat unit), which 
is a direct consequence of the self-localization effects in 
organic semiconductors and the concomitant positive de- 
pendence of the mobility (/Lt) of electrons and holes on 
concentration (n), 0, 13, Finally, the mobility /i(ri) 
exhibits a strong dependence on the doping mode, and 
a much stronger positive dependence is in effect when 
near-by compensating (and lattice-polarizing) counter- 
ions are present, as is the case in electrochemical and 
chemical doping, than when the electronic charge carri- 
ers are introduced in, e.g., a field-effect mode [4, 8]. 

In this Letter, we develop a model describing the elec- 
trochemical conversion of an organic semiconductor film 
from the undoped non-conducting state to the doped con- 
ducting state, and compare the results with recently ac- 
quired experimental data. It is demonstrated that the 
observed doping front in light-emitting electrochemical 
cell (LEG) devices can be described with a model based 
on a combination of a strongly concentration-dependent 
mobility for the electronic charge carriers and the Nernst 
equation in the doped conducting regions. 

A typical LEG consists of a solid-state active mate- 
rial, comprising an intimate blend of a fluorescent con- 
jugated polymer and an electrolyte, positioned between 
two electrodes. When a voltage equal to or larger than 
the band-gap potential of the conjugated polymer is ap- 
plied between the two electrodes {V > Eg/e), balanced 
charge injection (electrons at the cathode and holes at the 
anode) into the conjugated polymer is facilitated by the 
preceding migratory motion of the electrochemically in- 
ert ions and the corresponding formation of electric dou- 
ble layers at the two electrode interfaces. The injected 
electronic charge carriers are subsequently in effect elec- 
trostatically neutralized by the compensatory motion of 
ions in a process termed electrochemical doping; p-type 
doping (i.e., injection of holes and compensation by an- 
ions) takes place at the positive anode and n-type doping 
(injection of electron and compensation by cations) at 
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the negative cathode, and after a "turn-on time" a light- 
emitting p-n junction is formed in the inter-electrode gap 
0) The doping of a fluorescent conjugated polymer 
has a "dark" optical signature in that the formation of 
dopants (polarons) is concomitant with the formation of 
low-energy sites with limited radiative-decay probability 
from the excited state {i.e., quenched fluorescence). Ac- 
cordingly, by exposing an LEG device during operation 
in a dark environment to UV light (which excites the 
fluorescence of the conjugated polymer), it is possible to 
correlate the formation of a doped conjugated polymer 
to the emergence of dark regions with quenched fluo- 
rescence. Figure la presents a photograph of a planar 
LEG device during operation at F = 5 V (> Eg/e) and 
under UV light illumination in a dark room. The an- 
ode/cathode is marked with a "+/—". The fluorescence 
quenching of the "green" and the "red" components of 
the photograph along a representative path in the inter- 
electrode gap (as indicated by the dotted hne in Fig. lb) 
are presented in Figs. Ic and d, respectively. Two ob- 
served distinct dark regions with very low fluorescence 
intensity correspond to p-type and n-type doping. The 
regions originate at the anodic and cathodic interfaces 
and end with a very sharp and rather irregularly shaped 
front boundary. By correlating the size and growth of 
each doped region with the temporal evolution of the 
integrated current, it is possible to extract important 
information as regards to the doping concentration and 
profile. We have recently demonstrated that the doping 
concentration in the doped regions is very high, with a 
value of the order of ~10^^ dopants /m"^ (corresponding 
to '--^O.l dopants/repeat unit), [11|, and that the doping 
concentration behind the doping front is relatively con- 
stant [ilHQIii. 

The organic semiconductor material will be modeled 
as consisting of electrons and holes together with posi- 
tive and negative ions in a polymer structure. Diffusion 
and mobility of charges are of primary importance, which 
may be studied using the methods of plasma physics. 
Similar models have been presented before, but we will 
show that very important aspects have previously been 
overlooked. In particular, we include an injection barrier 
in the form of the Nernst potential in conjunction with a 
concentration dependent mobility, which are necessary to 
obtain agreement between theory and experiments. The 
particle species in the material satisfy the force balance 
equation 

nmv = — (7nV(0 — 4>n) — ksTVu — T^^nmw, (1) 

where n denotes the concentration, m is the mass, v is 
the velocity, q is the charge, ks is Boltzmann's constant, 
T is the temperature, r is the time between collisions, 
and 4> — (j)N is the overpotential experienced by the parti- 
cle. Here we have introduced a phenomenological injec- 
tion barrier (pN due to the difference in Fermi levels for 
the electrons and holes between the doped and undoped 




FIG. 1: (a) Photograph of a planar Au/{MEH- 
PPV+PEO-I-KCF3SO3/AU device with a 1 mm inter- 
electrode gap during operation at = 5 V and T = 360 
K. The anode and cathode are indicated by a {+) and a 
(-) sign, respectively. The device was operated in a dark 
room under UV illumination, so that the p-type and n-type 
doping regions are visualized as dark regions originating 
from the anodic and cathodic interfaces, respectively, (b) 
Enlarged section of the photograph in (a), as indicated by 
the solid line. The dotted line in (b) marks the path in the 
inter-electrode gap along which the fluorescence quenching 
data of (c) the green component and (d) the red component 
were extracted from the photograph. 

regions. Ions do not experience any barrier, so that (pN 
should be taken zero for ions. In the electrochemical sys- 
tem under study, the injection process is not limiting [iq . 
i.e., the electronic species on both sides of the barrier are 
at quasi-equilibrium, and ^at is chosen to be the Nernst 
potential 

(j>N{nh,e,T) = (fcsr/g/i,e)ln[7i/i,e/(n/i,e.oo -n/i,e)], (2) 

where nh^e,oo is the maximum concentration of holes (la- 
bel "/i") or electrons (label "e"). The Nernst potential 
originates from quantum effects, and the interpretation 
of in terms of an electrostatic force should therefore 
be done with caution. The Nernst potential is valid only 
in the highly doped region, and it should be set formally 
to zero in the undoped region, (f)iq{nh,e^T) = 0. When 
going over from one region to another, we match the so- 
lutions by fitting the functions and their derivatives. 

On the timescales of interest here, i.e., much longer 
than the collision time t, the collision force dominates 
over the particle acceleration mndv/dt, such that the 
left hand side of Eq. ([T]) can be neglected. Then v — 
—Dn^^'Vn ± i-LV{(j) — <^Ar), where the mobility is given 
by jj. — T \q\ /to, the diffusion coefficient D follows from 
the Einstein relation = \q\D/kBT and ± correspond 
to positive and negative charges. Substituting v into the 



continuity equations for all species we obtain the equa- 
tions for the semiconductor system 

dtrih.e - V • [±^h,en/i,eV((/) - 4>n) + Dh,e^nh,e] = 0,(3) 

dtn±-V ■[±ii±n±V(j) + D±Vn±] = Q. (4) 

The labels denote holes (ft.), electrons (e), positive (+) 
and negative (— ) ions. The electrostatic potential is 
determined by Poisson's equation £oV^(/) — —e{n^ + nh — 
n_ — He). Still, in our system, the quasi- neutrality condi- 
tion rt+ + nil — ri- — Tie = is satisfied with very good ac- 
curacy. We note that the quasi-neutrality condition does 
not imply a constant electric field strength, but rather 
indicates that small charge imbalance causes extremely 
large electric fields |17[. 

In what follows we will focus on the dynamics and 
structure of the p-type doping, but the analysis could 
of course equally well be applied to the n-type doping. 
Equations ([3]) and (|4]) admit a solution in a form of a 
localized planar doping front resembling a shock, moving 
with velocity U and converting the undoped polymer into 
a doped polymer (for a similar observation in ferroelectric 
crystals, see Ref. [3|). In the undoped region (label 0), 
we have a uniform electric field —dx4>Q, uniform initial 
concentration of ions n_|_o = n_o = no, and no holes, 
nho — 0. Behind the front, in the highly p-doped region 
(label cxd), the hole concentration is denoted by n/ioo, 
while the electric field —dx4>oa and ion concentrations 
n-too are obtained from Eq. ^ and the quasi-neutrality 
condition. In the reference frame of the doping front 
{dtn = 0), we may integrate Eqs. ^ and ^ to obtain 

DhOxUh = -UhU - nh^J,h^x{(t> ~ (t>N), (5) 
D±dxn± = {no - n±)U ± fi± {n±dx4> - mdx4>o) ,(6) 

for the holes and ions, respectively. Noting that the con- 
centration gradients vanish on each side of the doping 
front, we obtain the p-front velocity 
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U = -(no/n/ioo)(/-i+ + fi-)dx(f>o. 



(7) 



The electric field in the p-doped region is negligible due 
to the high mobility of the holes relative the ions. Diffu- 
sion does not influence the front velocity ([7]), but instead 
determines the structure and width of the front. 

A dimensional analysis suggests the characteristic 
width of the front Lf = D-/U. However, Lf does not 
portray the full structure of the doping front, since there 
are several different characteristic length scales within 
the front, from the undoped to the doped region. The 
hole mobility is highly sensitive to the concentration, and 
in the undoped region Hho/f^± ^ 1- Thus, from Eqs. ([5]) 
and ([6]) we also find that the characteristic length scale in 
the undoped region is given by (^ho//J.±)i/, which may 
be orders of magnitude smaller than Lf. As a result, 
one can expect a sharp head of the doping front with 
strong concentration gradients. On the other hand, for 
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FIG. 2: Model data for the concentration of ions as a func- 
tion of distance from the positive anode, resulting from the 
employment of: (a) concentration- depencieni mobility for the 
electronic charge carriers, (b) concentration- inciependeni mo- 
bility and the Nernst equation in the doped regions, and (c) 
concentration- dependent mobility and the Nernst equation in 
the doped region. 



the doped region, Eqs. (O, (HI) predict a slow variation of 
the front due to Hhoo/fJ'± ^ 1 and the Nernst potential. 
The concentration approaches the final value asymptoti- 
cally according to the power law (/i+ /fj,^ + l)Lf/{x — Xf ), 
where x — Xf measures the distance from the leading edge 
of the doping front. This gives a very smooth behavior 
of the doping front in the highly conducting part of the 
polymer. The analytical reasoning has been supported 
by the numerical solution to Eqs. (O, ([6]), described be- 
low and presented in Figs. 2, 3. Thus, the structure of the 
doping front and the generic features of the doping pro- 
cess are reproduced by our model, both analytically and 
numerically. Some of the length scales presented above 
could be measured in future high resolution experiments. 

As stated above, our model makes use of the concen- 
tration dependence of the hole mobility (see, e.g., Refs. 
0,i). In the numerical solution we employ the following 
empirical concentration-dependent mobility for holes as 
extracted from Ref. M 
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(8) 

Note that we have renormalized the mobility data in ac- 
cordance with data from Ref. to account for the lower 
mobility of MEH-PPV in comparison to the conjugated 
polymer under study in Ref. [20|. To first approxima- 
tion the ion mobility does not depend on concentration. 
To analyze Eqs. ^ an d (El) numerically, we use Eq. ([8|) 
and empirical data 12, Tj^llSl at the early stage of the 

d(f>o/dx = 3V/mm 
= 2.2 X 10-1° mVVs, 



doping process: uq = 3.1 x 10 



I.O X 10-1° mVVs 
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FIG. 3: Model data showing concentration of holes, and the 
overpotential profile as a function of distance from positive 
electrode. 



nhoo — 8-8 X iU'"ni and T = 360 K. The respective 
numerical solution is shown in Figs. 2 and 3. Figure 2 
shows the ionic concentrations as a function of distance 
from the anode (x) when the p-type doping front has 
progressed to x = 0.3 mm. We note the highly resolved 
X-axis. Our model solution is shown in Fig. 2c, while 
Figs. 2a and b illustrate the unphysical results with- 
out employing the Nernst equation and concentration- 
dependent mobility, respectively. First, by omitting the 
Nernst term, we obtain correct concentrations in the un- 
doped region but cannot achieve physically sound results 
in the doped region (Fig. 2a) . Second, if we instead keep 
the Nernst term, and assume a constant hole-mobility, 
the doped side of the front is described quite well. How- 
ever, ionic concentration levels diverge rapidly in the un- 
doped region (Fig. 2b). If we, on the other hand, in- 
clude both the correct concentration dependent mobility 
(IH) as well as the Nernst term ([2]), we obtain the physi- 
cally sound result shown in Fig. 2c. Here we can see the 
emergence of a doping front that has a shape in good 
agreement with the experimental observations of Fig. 1 . 
The doping front in Fig. 2c has a characteristic width 
~ 0.1 — 0.2 mm, a consequence of the long tail due to the 
Nernst term in Eqs. ^ and ([7]). At the leading edge of 
the front, the gradients look very sharp as compared to 
the trailing edge. This follows from the low hole mobil- 
ity in the undoped region, and a higher resolution of the 
leading edge demonstrates that the characteristic length 
scale is {fj,ho/n±)Lf w 2 x 10~^mm, see the insert of Fig. 
3. Figure 3 shows the concentration of holes within the 
front, and the overpotential profile as a function of dis- 
tance from the anode. Again, in agreement with the pre- 
vious discussion, we observe a very sharp leading edge of 
the front and a smooth long tail. The plot demonstrates 
a weakly varying over-potential at the trailing edge of 
the front. We also compare quantitatively the theoretical 
predictions for the doping front velocity and the experi- 
mental measurements. Using Eq. ^ and the experimen- 



tal data we calculate the front velocity at the beginning 
of the doping process as U = 13.5 x 10~^m/s. The ex- 
perimental measurements provided the same initial value 
for the front velocity, U = 13.5 x 10~^m/s, within the 
accuracy of measurements. As the doping progresses, the 
p- and n-type doping fronts accelerate toward each other 
|la |. due to the decrease in distance between the fronts. 

Thus, our model is in very good agreement with the 
experiments, both qualitatively and quantitatively. It 
predicts the quantitative value of the doping front ve- 
locity with good accuracy. The model also gives good 
agreement between the observed and calculated front 
thickness. Furthermore, it reproduces the correct qual- 
itative front structure as a combination of smooth tail 
and very sharp leading edge. Though the present anal- 
ysis is one-dimensional, it can be easily extended to a 
multi-dimensional case. However, the study of multi- 
dimensional effects is left for future research. 
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